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In this paper Maxwell's partial differential equations and the boundary 
conditions for waveguides fdled with a heterogeneous and non-isotropic 
medium are converted into an infinite system of ordinary differential equa- 
tions. This system represents a generalization of "telegraphist's equations" 
for a single mode transmission to the case of multiple mode transmission. 
A similar set of equations is obtained for spherical waves. Although such 
generalized telegraphist's equations are very complicated, it is very likely 
that useful residts can be obtained by an appropriate modal analysis. 

From a purely mathematical point of view the problem of electro- 
magnetic wave propagation inside a metal waveguide reduces to obtain- 
ing that solution of Maxwell's equations which satisfies certain boundary 
conditions along the waveguide and certain terminal conditions at the 
ends of the waveguide. If the medium inside the waveguide is homo- 
geneous and isotropic and if the cross-section of the waveguide is either 
rectangular or circular or elliptic, the desired solution is obtained by the 
method of separating the variables. The method works for some other 
special cross-sections. It works also if the medium inside a rectangular 
waveguide consists of homogeneous, isotropic strata parallel to one of 
its faces. Similarly, it works if the medium inside a circular waveguide 
consists of coaxial, homogeneous, isotropic layers. But in general if the 
medium is either nonhomogeneous or non-isotropic or both, the method 
does not work. The mathematical reason for this is that the solution is 
of a more complicated form than a simple production of functions, each 
depending on a single coordinate. Any function that one usually en- 
counters in physical problems, and therefore a solution of Maxwell's 
(Miuations, may lie expanded in a scries of orthogonal functions. Sets of 
such functions are provided by the solutions for waveguides filled with 
homogeneous media. Such functions already satisfy the proper boundary 
conditions and the problem is to obtain series which also satisfy 
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Maxwell's equations. From the physical point of view this method 
represents a conversion of Maxwell's equations into generalized "tele- 
graphist's equations." 

Thus it is already known that Maxwell's partial differential equations 
and the boundary conditions along a waveguide are convertible into a 
set of independent ordinary differential equations, each resembling tele- 
graphist's equations for electric transmission lines. 1 Each equation de- 
scribes a "mode of propagation" in terms of concepts well known in 
electric circuit theory. A waveguide can be considered as an infinite 
system of transmission lines. If the medium inside the waveguide is 
homogeneous and isotropic and if the surface impedance of the boundary 
is zero, the method of separating the variables enables us to obtain a set 
of "normal", that is, uncoupled modes of propagation. Any irregularity 
or "discontinuity" in the waveguide provides a coupling between 
some, or all, modes of propagation. The irregularity may be in a boundary 
of the waveguide or in the dielectric within it. A heterogeneous dielectric 
may be considered as a homogeneous dielectric with distributed irregu- 
larities. 2 Similarly a heterogeneous non-isotropic dielectric may be con- 
sidered as a homogeneous isotropic dielectric with distributed irregu- 
larities. Such irregularities provide a distributed coupling between the 
various modes appropriate to homogeneous isotropic waveguides. Our 
problem is to calculate the coupling coefficients. The generalized tele- 
graphist's equations, obtained in this manner, are very complicated in 
that they represent an infinite number of coupled transmission modes. 
They are useful, however, in suggesting a physical picture of wave 
propagation under complicated conditions, and can be used in approxi- 
mate analysis when we can ignore all but the most tightly coupled 
modes. For example, this picture was successfully employed by Alber- 
sheim 3 in studying the effect of gentle bending of a waveguide on propa- 
gation of circular electric waves. In this caseit was important to consider 
the coupling between only two modes, TE i and TMn , which have the 
same cutoff frequency in a straight waveguide. More recently, Stevenson 
obtained exact equations for waves in horns of arbitrary shape. 4 His 
equations express the propagation of the axial components of E and H. 
The various modes are coupled through the boundary of the horn. In 

1 S. A. Schelkunoff, "Transmission Theory of Plane Electromagnetic Waves," 
Proc. Inst. Radio Engrs., Nov. 1937, pp. 1457-1492. 

2 S. A. Schelkunoff, "Electromagnetic Waves," D. van Nostrand Co., (1943), 
pp. 92-93. 

3 W. J. Albersheim, "Propagation of TE i Waves in Curved Waveguides," 
Bell System Tech. J., Jan. 1949, pp. 1-32. 

4 A. F. Stevenson, "CJeneral Theory of Electromagnetic Horns," J. Appl. 
Phys., Dec. 1951, pp. 1447-1460. 
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the present paper we shall consider waveguides of constant cross-section 
and conical horns of arbitrary shape rilled Avith a heterogeneous and 
non-isotropic dielectric and derive the equations for propagation of the 
generalized voltages and currents representing the transverse field com- 
ponents. The various modes are coupled through the medium. It is very 
likely that our equations can be generalized to include the coupling 
through the boundary. 

To understand the mechanism of coupling between the various modes 
through the medium consider Maxwell's equations 

curl E = -jaB, curl H = "J + j<»D, (1) 

where e J is the density of conduction current while the other letter 
symbols have the usual meanings. In the most general linear case the 
components of B and D are linear functions of the components of H 
and E respectively, with the coefficients depending on the coordinates. 
These equations can always be rewritten as follows 

curl E = -jowH - M , curl H = jcoeE + J, (2) 

where M and / are the densities of magnetic and electric polarization 
currents. 

M = MB - uH), J = C J + ju(D - eE), (3) 

and n, e are constants (not necessarily those of vacuum). If M and J 
were given, they would act as sources exciting various modes of propa- 
gation in a homogeneous, isotropic waveguide. If M and J are functions 
of H and E, they can still be considered as the sources, acting on power 
borrowed from the wave, of the various modes. Thus M and J will 
provide the coupling between the modes existing in a homogeneous, 
isotropic waveguide. 

Thus in order to derive the generalized telegraphist's equations we 
shall first consider the various modes of propagation in a homogeneous 
isotropic wave guide. Each mode is described by a transverse field distri- 
bution pattern 6 T(u, v), where u and v are orthogonal coordinates of a 
point in a typical cross-section. This function is a solution of the follow- 
ing two-dimensional partial differential equation 



AT = 



!f» (*<*)+■• (5 ^Yl.-vr, (4) 

l e-» [_du \ei du/ dv \e 2 dv / J 



5 Sec Reference 2. 

6 S. A. Schelkunoff, "Electromagnetic Waves," D. van Nostrand Co. (1943), 
Chapter 10. 
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where x is the separation constant and e r , e 2 are the scale factors in 
the expression for the elementary distance 

ds = ei du + el dv~. (5) 

In the case of TM waves the T-function must vanish on the boundary 
of zero impedance. This boundary condition restricts x to a doubly in- 
finite set of values Xmn with the corresponding functions T mn . In the 
case of TE waves the normal derivative of the ^-function must vanish 
on the boundary of zero impedance. Since we have to consider both 
types of waves simultaneously, we shall distinguish between them by 
enclosing the subscripts in parentheses for TM waves and in brackets 
for TE waves. The double subscript designation of various modes has 
been standardized only for rectangular and circular waveguides. For 
waveguides of other shapes the standard is to use a single subscript by 
arranging the modes in the order of their cutoff frequencies. For con- 
venience, we shall use this convention in the general case and denote 
TM modes by T( n )(u, v), and TE modes by r l\ n ](u, v). The correspond- 
ing cutoff constants will be xw and xm • In what follows it is under- 
stood that whenever the T-functions should be designated by double 
subscripts, our single letter subscripts should be considered as symbols 
for ordered double subscripts. 

The transverse field components may be derived from the potential 
and stream functions, V and II for TM waves and U and ^ for TE 
waves. Thus 

E t = - grad V - flux *, H t = flux II - grad U, (G) 

where the components of the gradient and flux of a scalar function W 
are 

grad u W = — — , grad, W = — — , 

e x ou e 2 ov 

(7) 

fl W dW fl TIT" dW 

flux,, W = , flux,, W = — — — . 

e> dv ei ou 

The T-i unctions corresponding to the various modes of the same va- 
riety are orthogonal ; that is, the following integrals over the cross-section 
vanish, 

JJT (n) T (m) dS = 0, fJT M T lm] dS = 0, tim* n. (8) 

It should be stressed that T( n) and T\ m \ are not, in general, orthogonal. 
7 See Reference 6. 



(9) 
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Similarly the gradients of the ^-functions of the same variety as well as 
the fluxes, are orthogonal, 

//(grad r ( „,)-(grad T ( „,) dS - //(flux T (n) )-(flux T (m) ) dS 

- //(grad r w ) • (grad T lm] ) dS = //(flux I» w ) • (flux T lm] ) dS = 0, 

if m ^ n. The following gradients and fluxes of the ^-functions are ortho- 
gonal for all m and n, 

/"((grad Too) -(flux T [n] ) dS = //(grad T [n] )-(nux F (w) ) dS 
JJ (10) 

= //(grad T (n) )-(flux T (m) ) dS = 0. 

On the other hand, grad T [m] and flux T M are not, in general, orthogonal. 
If all modes are present, the potential and stream functions are 

V = -V in) (z)T (n) (u,v), n = -i w (*)r©o(«,tO, 

* = -7 w (2)r w (w, v), u = -J w (*)r w («, »), 

where the tensor summation convention is used: whenever the same 
letter subscript is used in a product, it should receive all values in a 
given set and the resulting products should be added. The negative 
signs have been inserted in order to avoid a preponderance of negative 
signs in later equations. Substituting in (9), we have 

E t = 7 (n) grad T in) + V [n} flux T M , 

(12) 
H t = -I M flux r< w) + I ln] grad T [n] . 

The !F-functions for the various modes are determined by equation 
(4) and the boundary conditions except for arbitrary factors related to 
the power levels of the modes. If we choose these constants in such a 
way that 

//(grad T) ■ (grad T) dS = x 2 JJt* dS = 1, (13) 

then the complex power carried by the wave is given by an expression 
similar to that in an ordinary transmission line, 
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Hence, the V's and 7's correspond to the voltages and currents in coupled 
transmission lines. 

In an expanded form equations (12) are 



—= 1 V[n] — , £iv = K(„) — — - V [nj — — , 

ei du e-> dv e 2 dv e\ au 



TJ J dT(„) j dT [n ] IT — T dT(„) j dT[n\ 

To these we add the following expansions for the longitudinal compo- 
nents of E and H 

E z = x M V t , M (z)T M (u, v), H z = XMliMMTwiu, v). (16) 

Equations of this form satisfy automatically the boundary conditions 
on Ei and H z . The multipliers x« have been inserted arbitrarily in order 
to make the physical dimensions of the second factors to correspond to 
those of voltage and current. 

Let us now write Maxwell's equations in an expanded form 

dE, dE v . D dH, dH„ . n 

= — 70)/),, , — — — = ]<j)JJ u . 

e 2 dv dz J e-idv dz J 

dE u dE, . dH H dH z > 

-r- - — -r- = —jutt. , -T- - — r- = JuVv , (17; 

dz e x du dz ei du 

d(e 2 E v ) d(eJ<J u ) _ d(e 2 H r ) d(eiH u ) n 

du dv J ' du dv J 

Substituting from (15) and (lfi) in the left column of (17), we find 

dT M dV (n) dT w dV ln] dT [n] 

X(n) V:.(n) — ~ — "J — + —J "— T— = —jwB,, , (18) 

c 2 dv dz Co dv dz ei du 

T/ dT M dV M dT w . dV [n] dT ln] . _ . . 

<rr (n) _ v d /c2 d'l^A _ v aX, _ F d /e, ar [n ,\ 

K(n) du dv [n] du U du )~ (n) dv du [n] dv \e» dv J 

(20) 

In view of (4) the last equation reduces to 

xl\V lnl T [n] = -j*B z . (21) 
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Multiplying (18) by [-dT {m) /e 2 dv] dS, (19) by [dT^/c, du] dS, adding, 
and integrating over the cross-section, we obtain 

_ X(m)7 „ ( „, + <T« = iw // ( B . »jfe - B. g&) « . (22) 

' In the first term the summation convention should be ignored. Multiply- 
ing (18) by [8T lm] /ei du] dS, (19) by [dT lm] /e 2 dv] dS, adding, and in- 
tegrating we find 

•Si = -j a ff ( Bm 3 ^! + B. d -%i) dS. (23) 

dz J J \ ej du e 2 dv / 

Multiplying (21) by T [m] dS and integrating, we have 

V [m] = -jm JJb.T w dS. (24) 

Subjecting the right column of (17) to a similar treatment, we obtain 
three additional equations. Summarizing, we have 

'-Zp = i« ff (B u ^ - B. ??W) dS + xw^w . (25) 
dz J J \ e 2 dv e y du J 

d Jp = - ju ff ( Dh %> + D. *?W) dS, (26) 

dz JJ \ Ci du e 2 dv / 

d -Y^l = -jco ff (b u ^ + B, a %l) dS, (27) 

dz J J \ e\ du e> dv / 

S=i = jw ff f- Du W=* + D, ??W) dS + XM/..M , (28) 

dz J J \ e>dv e, du/ 

V lm] = -j<* fJB t T lm] dS, I (m) = -jw fJD a T im) dS. (29) 

In the last terms of equations (25) and (28) the summation convention 
should be ignored. 

If the components of B and D are linear functions of the components 
of H and E respectively, then with the aid of (15) and (16) they can be 
expressed as linear functions -of V< n ) , F(»j , /<„> , /[„] , V m ,m , /*,[„] - 
Solving (29) for V : , M and I z , [n ] and making the appropriate substitu- 
tions in (25), (26), (27), (28), we obtain the generalized telegraphist's 
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equations in the following form 

d'(m) rr T V T r T ~\T Vr F 1/ 

z = — Z(m)(n)i(n) ~ ^(m)[n]l[n] ~ I (m)(n) » (n) ~ *■ (m)\n\y [n] , 

— = — 1 (m)(n) » (It) ~ ' (»)W » ["1 _ J fm)(n)i(n) i (m)[n]i[n] , 

d " (30) 

r^ = — ^[m][n]/[n] ~ Z[ m ] („)/(„) — T[ m ] [„] !"[»] — T[ m ]( n )V {n) , 

az 

dl 



f^ = —Y [m ][n]V [n ) — Y [m ) M V M — T [m ][ n )I[n) — T[ m }(n)Iin). 



dz 

The transfer impedances Z, the transfer admittances Y, the voltage 
transfer coefficients V T, and the current transfer coefficients T T between 
various modes are in general functions of z. They are constants if the 
properties of the waveguide are independent of the distance along it; 
in this case the problem of solving the generalized telegraphist's equa- 
tions reduces to solving an infinite system of linear algebraic equations 
and the corresponding characteristic equation. 

Similar equations may be derived for spherical- waves either in an un- 
limited medium or in a medium bounded by a perfectly conducting coni- 
cal surface of arbitrary cross-section. If the latter is circular and if the 
flare angle is 180°, we have a plane boundary. Hence, the case of spheri- 
cal waves in a non-homogeneous medium is included. In the spherical 
case we shall use the general orthogonal system of coordinates (r, u, v) 
where r is the distance from the center and (u, v) are orthogonal angular 
coordinates. In this system the elements of length ds and area dS are 
given by 
(l / = d / + r \ e \ f / (i 2 + e l dif), dS = AM, dn = 44 du do. (31) 

The transverse field components may be expressed in a form similar 
to that for waveguides 

rE t = - grad V - flux II, rH, = flux II - grad U, (32) 

where grad and flux of a typical scalar function are defined by equations 
(10). Instead of (11) we have 

V = -V M (r)'l\ n) (u, >'), n = -I M (r)T w (% v), 

* = -V M (r)T M (u,v), U= -I [n] (r)T [n] (u,v), 

where the T-functions satisfy equation (4) and appropriate boundary 
conditions. These functions, their gradients and fluxes are orthogonal. 
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They are assumed to be normalized as follows 

jf (grad T) • (grad T) dSl = x 2 JJt 2 dQ = 1 , (34) 

where dJ2 is an elementary solid angle. Hence, equation (14) will again 
represent the complex power flow in the direction of propagation. 
The various field components may then be expressed as follows 

r& u = K ( „) — h K [n] — — , r£,„ = K< n) — — — V [ni — — , 

e x du e 2 dv e 2 dv c t du 

~n — t d7V) , T dT [n] „ , dT {n) T dT [n] (35) 

rH„ = /(„) — — - + I[n) — — , rti u = — /(„) — — + i [n] — — -, 
ei 5m e 2 dv e 2 3u e t du 

T E r = X(n) V r ,(n) ^(n) , T H T = X[n] /f.Jnl T[„j . 

It should be noted that the physical dimensions of V T ,(„) and I r ,in] are 
not those of voltage and current. Substituting in Maxwell's equations 
and using transformations similar to those in the case of plane waves, 
we find 



dV 



(HI) 



dr 



dl(m) 

dV [m ] 
dr 



J J \ e\ au e-> dv / 

*. /if (rB. ?5s + rB. £&) dO, (36) 

J J \ ei du e 2 dv / 



dI M _ „• ff ( ' „n dT w _l -n dT i"*! bo±v. , v ~h 



, =jco(- rD u ^ + rD e ^p dn + x[.ir-Vr iW , 
dr J J \ e 2 dv ej 3m / 

V lm] = -jco ff(r 2 B r )T lm] dn, /(«, = -> ff(r 2 Dr)T {m) (m. 

Returning to the plane wave case and assuming the following general 
linear relations 

B u = UuuH u + fi uv H v + IJ-uzHz , D u = e uu E u + Cut^ti + e uz S r , 

£„ = M.»ff« + HwH v + mJZ, , Z>„ = *JF« + *J? f + e vt E z , (37) 

Bz = u, u H u + ft-ff- + m«# 2 , D z = « 2u tf u + e, v E v + «J0. , 
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B u = I 



B 



= Art -I 



we find 

R .T |"-« dT <">4-u dT «l 4- 7r , [a ^ + U ^l 

dT w dT (n) ~\ [ dT M dT M l 

-^ + M " ^ J + [nl L M -' U ^ + M " ^ J 

+ 7z.[nJ/izzX[n]T' [n J , (38) 

L Ci5m e2avj L Biol) Ci 3u_ 

+ Vz,( n )€ uz X(n)T( n ) , 

D » = F(n) L e,u ?aS + *• ^7a7 J + 7[nl _ € " tt -£Fv ~ €vv i^ J 

+ V z,{n)*vzX{n)T (n) , 

[_ 0i aw e2 5v J L ^2 dv ei du J 

+ V z ,( n )€ Z zX(n)T( n ) . 

Substituting from equations (38) into equations (25) to (29) we ob- 
tain 

dV ™ it [f[ dT ™ ar « 4- „ dT ™ dT ™ 

IT = "• 2w7( " ) // L Muu ^ ^aT + ^ ^au 1^5 

d dw e 2 av e 2 av ei du J 

+ jo}i [n] jj \n uu — z — ' — — - — - v' 



e\ du e 2 dv 



ei dv ei du 



dF[ n ] dT( m ) dT[ n ] dT( m )\ ,„ 

e2 dv e 2 dv ei du e\ du J 
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dlt 



llfm) • XT- ff 

W = - JuV ™ Si 



dT w dT (m) dT w dT (m) 



e\ du d du 



e 2 dv <?2 dv 



dT( n ) dT( m ) , ' dT(„) dT( m ) 

T Cm- ^ ~ r c < 



e 2 dy Ci du 



ei du e 2 dv 



dS 



+ *F W // 



- *F..„ // [*. 






e 2 dv ei dw 



ei dw e 2 dv 



+ 6„, 



ei 5m 



+ 6, 



32'[ n ] dT( m ) 
dT 



BT M dT 



(m) 



e 2 dv Ca dv 



(40) 



r/.S' 



e 2 



8^-^ 



d7 



(fe 



M-irf« L uu 



M.t 



djT(„) d!T[ m ] 



6T(„) d ?'[„,] 

e 2 dv ei 3m *"" ei dw e 2 dy 

3T W dT [m] 



e\ du C\ du 



+ 



- j<*I M // 



djVl 3r [W ] . djT [n ] dT[m] 



e\du ei dit 



e 2 dw e 2 9t> 



di(„) dT[„,] ] , Q 

Mru — r— do 

e 2 dy e 2 dv J 

(41) 



• ay J 



e 2 dv e! dw 



ei dw e 2 



Ihu 



dT 



['»] 



— e UM 



d dw 

d2\„) dT^,,,] 



e 2 dv J 

3T(„) dT [m j 



+ e„„ 



d dw e 2 5?; e 2 dv ei du 

dT( n ) dT[ m ] 



^uv 



- **« // [ 



e 2 dv e 2 dv 



+ e 



dT[„] 3T[m] ^^[r.] dT[ m ] 

€i,n : : — ~\~ Gvr 



e 2 dv e 2 dy 



ei dw ei dw 



32jw] 9T[w] _ dT [n ] dT[ m ] 



d dw e 2 dv 



d^(n) 37jm] 1 j„ 

"" ei dw ei dw J 

(42) 

dS 



e 2 dv e\ du _ 



+ i«7,. (M) ff I"-*. ^% ] + e M ^1 XM T W dS + X[»]/..t«] , 

J J |_ e 2 ov ei aw J 

^.[nl /J J0>HzzXln]T [n] T [m] dS = —V [m ] 
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(»>> 



ei du e-> dv 



T im) dS (44) 



- v [p) II Ju\ — e*« — — + t 



Co di- 



et du 



T( m ) dS. 



If we solve the last two equations for I z ,[ n ] and F r ,(„) and substitute 
in the preceding four equations, we shall obtain the telegraphist's equa- 
tions in their final form (30). Thus, let 



(45) 



Z[mHn] = 1 1 J0>HzzXln]T [n ]T [m l dS, 



r F ( „ (, l ) = jj J0)€ zz X(n)T in) 7\ m) dS. 

From these coefficients we obtain another set 

: Z[„][ m ] = normalized co-factor of z Z[ m ][„] 
r Z(n)(»o = normalized co-factor of Z F< „,)(„> 

Then, 

Iz,[»] = —V [»>]'¥[■>][»>) 



(46) 



d7 (p ) dl( P ) - „. 

r- — Mz» — r- i [m] do 
2 3?' ei dif 



ff ; T d T [v] dT[p]\ 7 , ,., 

J J Ci dw ei dv J 



V z ,(n) = — I \m)"Z („)(»i) 

— V(p{Z ( n) ( m ) 

+ V[p] Z Z(n)(m) 



(47) 



;« 



jw 



d?0») I d7\„) 



e 2 dv 



eo dv 

dT [p ] 
ei du 



T(,n) dS 



T(m) dS. 



Before substituting in equations (39) to (42), the summation index m 
in (47) should be changed to avoid conflict with in in the former equa- 
tions. It does not seem necessary to make these final substitutions in 
their most general form. The results are very complicated and in prac- 
tice the various coefficients are not independent. Some coefficients may 
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vanish; others may be small. In isotropic media, Mu« = /*»» — M** = 
n, e uu = t vv = e lz = e and the mutual coefficients vanish. In gyromag- 
netic media subjected to a strong magnetic field in the z-direction, the 
permeability coefficients of superposed ac fields are 

M«« = M«w = M> M»u — —Puv, M*u = Hzv = Vuz = /*«* — 0. (48) 

If the entire waveguide is filled with such a medium, assumed to be 
homogeneous, equations (43) and (44) become 



I*,M i W M«X[»] // T[ n ]T[m]dS = — V[ m ) , 

Vz.wJueXin) J J T( n) T im) dS = —/(*») . 



(49) 



In view of the orthogonality of the T-functions and the normalization 
conditions (13), we have 

I..M = -£=*- V lm] , V 2 , (m > = -^ / (m) , (50) 

where the summation convention is waived. In this case all the transfer 
coefficients in equations (30) vanish, 

V T( m )(„) = V T(_ m )[n] = T [m ][„] = T[ m ](„) = T{ m )(n) = T( m )[n] (51) 

= T[m]M = T [m ][ n ) = 0. 

The transfer impedances and admittances are 
£(«,)(„) =0, if n ^ m, 

2 

= jojju + 4^ } , if n = m; 
jcoe 

ff [dT [n] dT, m) . dT ln] dT (m) ~] , , 
Zc«^t«i = — 7wu„„ // — — — — - H — — 7-* le&dudv; 

7(m)(») = 0, if w 5^ m, 

= jue, if n = m; 
F WW = 0, all m, n; (52) 

= jw/i, if w = ra; 



' C. L. Hogan, "The Ferromagnetic Faraday Effect at Microwave Frequencies 
and Its Applications — The Microwave Gyrator, Bell System Tech. J., Jan. 1952, 
p. 9. 
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_ ff\ dT( n ) dT[ m ] . dT( n ) dT[ m ] - , 

Zmm = jwuv // — -— -—- + — t 1 — \r wi du dv; 

J J \_e } du ei du e 2 dv e 2 dv J 
Y[m][n] = 0, if n ?± m, 

2 

= ^ + M if n = m; 

Fhh = 0, all m, n. 

We note that Z (m)[ni = -Z ln] (_, n) ; Z [m][n] = -Z [n][m] , (n ^m). 

In rectangular waveguides we choose cartesian coordinates; 
then C\ = Co = 1, u = x, v = y and 

- x . , N -l/2 • Ptx . qiry 
T( Pq) = 1 P3 X( P \)W sin H — sin ±f- , 

a b 

_x . , N -i/2 strx try ,__>. 

fwi = UxMab) cos — cos -^ , (53) 

2 2 2 2 

2 _ 2 _ p t , qv _ 2 

X(PQ) — X[P<?] °- ~^ I , 2 — Xpg > 

where l pg = 2 if neither p nor q is equal to zero and lo« = lpo = y/% 
Hence 

v> /"/"T/i./ > • s,ra; P"* ^ • w 

X // {b/a)sp sm — cos - — cos -r 2 sin 2— 
JJ a a b b 

dx dy 

W-wKt/a) 1 + ft/WTU - (-) 8+p ][i - (-) ?+ '] 

(54) 



-f- {a/b)tq cos — sm - — sin ~ cos ^ 
a a o o 



Ww x P ,x.U 2 -p 2 )(g 2 - « 2 ) 



if s t* p, q ?* t, 

= 0, if s = p or g = £; 

7 , ip,i.»(pV - gV)[i - (-r +p Ki - (-r + i 

$8 9* p,q 9* t, 

= 0, if s = p or g = t, but not if s = p and q = t, 
= jw/i, if s = p and 5 = t. 
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Some of the mutual impedances vanish; thus 

Z (pg)[ sn = 0, (55) 

if cither p + s or q + t is even. If p + s, as well as q + /, is odd, 

_ ■4-lp,l.jiw/i. y pg[(s/a) 2 + («/6) 2 ] / Kfi v 

XpaX,/(s 2 - p 2 )(? 2 - t-) 

Similarly 

%[pq][st] = 0) (57) 

?/ tt'tAer p + s or ry + £ is even, provided p ^ s and 9 3^ /. If p + s, as 
well as </ + /, is odd, 

4-lpgl.lCpV - qS 2 )jutl X y ,-qN 

1 XpvXstis 1 - pW - t-)ab 

Consider now the set of modes which includes TE [10 j. This set in- 
cludes TE [01 ] modes and all the other modes which are coupled to either 
of these modes. Noting that there are no TM (p0) and TM (0 ,) modes, we 
obtain the following table in which those modes which do not belong 
to the set are marked with a bar: 

TE[io] , TE[ou , 

TE [2 o] , TE[ii] , TE[02] , TMcu) , . . 

(59) 

TE [301 , TE [!tt ] , TE[i 2] , TE[os] , TM m , TM (12) , 



TE[4„] , TE [8 i] , TE |22] , TEtu] , TE [0 4] , TM (M) , TM (22) , TM(ii) , 
From the preceding equations we obtain the coupling impedances, 

o Q 

£[10] [01] = —juHxy , Z[01U10] = —-7,JUVxy, 

2(10] [30] = £[30] [10] = £[01] [03] = £[03] [01] = 0, 

g 
£[30] [01] = —Z[01] [30] = Z[10] [03] = —Z[03) [10] = ~~ o J^Hxy , 

£[2i] [io] = -£[io] [ 2 i] = -s-s-icWwU + 4(o/a) 2 ]" 1/2 , 

07T" 

£[oii[i2] = -£[i2, [ oi]'= ^JWx U [l + 4(a/6) 2 f 1/2 , (60) 

07T- 

£[10] [12] = £[1 2 ][10] = £[01] [21] = £[21] [01] ■= 0, 

1 I* ^O 

£(10]( 2 1) = — £(21)[10] = o . — JUPxyl* ~T \fl/U) J > 

07T" 
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Z [0 i](i2) = — £as>[oi] = Q s- i»mJ4 + [b/a)*T m , 

2(12) [10] = Z[10] (12) = ^(21) [01] = 2[01](21) = 0. 

The principal effect of the gyr< (magnetic medium on the TEuoj and TE [0 i] 
modes may be understood by taking into account their mutual coupling 
but ignoring their coupling to other modes. The equations of propaga- 
tion become 

, = —jo/if [io] — join X y(8/ir 2 )I[oi] , 
az 



m = -i*»+?-^)Vm 



dz \ joifizzd 

dV 

dz 

dl 



(61) 



[ — = jo)Mx tf (S/7r 2 )/[io] — jo}fil[oi] , 



(02) 



j~ = - [jo* + ~. rr, ) 7[oi]. 

i-Z \ JUVzzO 1 / 

For exponentially propagated waves we have 

V m = fiM]«"*, F [0 i] = t m e~ ifi \ 

•/[in] — l[iO]C , -/[Oi] — -MoiiC 

When the mutual permeability is zero, we have two independent modes 
whose phase constants are 

(2\l/2 / 2\l/2 

fl* ~ ^) , 0oi = (a, 2 /* " ^J • (63) 

The phase constants of the perturbed modes may be expressed in terms 
of the unperturbed constants and the coefficient of coupling. When the 
losses are neglected, the mutual permeability is a pure imaginary. In 
this case it is convenient to define a real coupling coefficient 

k = $%* . (64) 

Substituting from (62) in (61) and using (64), we find 



0V poi = wM?[ioj — jo}nkl[oi) , /3/[io] = I we mio] , 

\ (aHuOr/ 

( tT \ - 

/3F[oi] = jufMkl [io] + wm/[oi] , PI [oi] = ( toe — pil^ioi] 



(65) 
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Eliminating F[io] and V[ i\ , we have 

(0 2 - Pio)I[w] = — jkfivUm , 

(66) 

(/3 2 - j8oi)f [oi] = jkp&im ■ 

Multiplying term by term, we obtain the characteristic equation 

P A - (tfo + Pli)P 2 + (1 - fcVX = 0. (67) 

Solving, we have 

2 = iGA + /&) ± H(A -^oi) 2 + «WA] Vi . (68) 

The effect of coupling is to increase the larger phase constant and de- 
crease the smaller one ; that is, to make the slower wave slower, and the 
faster wave faster. 

Let us assume a > b; then ft > jSoi • Taking the upper sign in (68) 
and substituting in the second equation of the set (66), we have 

2[oi] j/cQ3oi//3io) 1 /fto _ 0oi\ , fiq x 

T^oi " V + (P 2 + &) W ' "2 W /3w/ ' * ' 

From (65) and (69) we find 

V"[on _ 0io ^lon jkifJio/pm) /„ v 

f[io] "" 2 i/[io] "" P + (P 2 + W* * 

Hence, the ratio of the power carried in the TE [0 u mode to that in the 
TE[ioj mode is 

Po\ V[0l)I[01] __ * (ni\ 

pTo " y[io]/fio, " [ P + (p 2 + fc 8 ) 1 ' 2 ]' " 

This ratio increases as k increases and p decreases. 

If the phase constants of the uncoupled modes are equal, then p = 
and Poi = Pan for all values of the coupling coefficient. In this case (68) 
becomes 

/3 2 = A(L ± k) or p = ft (l ± k) m . (72) 

In terms of the original constants, 

0=[(, ± l^)(^--^)f. (73) 

The cutoff frequencies of both normal modes are seen to be independent 
of either the transverse permeability or the mutual permeability. Since 



GENERALIZED TELEGRAPHIST'S EQUATIONS 801 

Hxv is a pure imaginary, it effectively increases the transverse permeabil- 
ity for one mode and decreases it for the other. 

To evaluate the effect of higher order TE and TM modes on wave 
propagation we may substitute from (68) in all terms of the character- 
istic equation for telegraphist's equations except the first two diagonal 
terms and recalculate the /3's. Alternatively we may replace TE[i ] and 
TE[o« modes by the normal modes just obtained, recalculate the cou- 
pling coefficients, and evaluate the effect of the mode with the greatest 
coupling to the modes under consideration. 



